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ABSTRACT 

We investigated the spectral evolution in the 2-700 keV energy band of Gamma-Ray Bursts 
(GRBs) detected by the Gamma-Ray Burst Monitor (GRBM) and localized with the Wide 
Field Cameras (WFCs) aboard the BeppoSAX satellite before May 1998. Most of them have 
been followed-up with the Narrow Field Instruments aboard the same satellite. In the light 
of these results we discuss open issues on the GRB phenomenon. We find that the optically 
thin synchrotron shock model (SSM) provides an acceptable representation of most of the time- 
resolved GRB spectra extending down to 2 keV, except in the initial phases of several bursts 
and during the whole duration of the quite strong GRB970111, where a low-energy photon 
depletion with respect to the thin SSM spectrum is observed. A strong and time variable low 
energy cut-off, consistent with absorption effect, is observed during the prompt emission of 
GRB980329. We find that the X-ray afterglow starts at about 50% of the GRB duration, and 
that its fluence, as computed from the WFC light curve, is consistent with the decay law found 
from the afterglow NFI observations. We also investigate the hydrodynamical evolution of the 
GRB in our sample and their associated afterglow, when it was detected. We find that the 
photon index of the latest spectrum of the GRB prompt emission is correlated with the index 
of the afterglow fading law, when available, as expected on the basis of an external shock of a 
relativistic fireball. We also find that for most of the GRBs in our sample the late emission is 
consistent with a slow cooling of the shock. Adiabatic shocks appear more likely than radiative 
shocks. Parameters of the shocks at earliest times have been derived. 

Subject headings: gamma rays: bursts — gamma rays: observations — X-rays: general — 
hydrodynamics: — shock waves 



1. INTRODUCTION 

In the last two years a big step forward has been 
accomplished in the GRB astronomy, thanks to the 



BeppoSAX( [Boella et al. 1997aD capability of precisely 
positioning these events soon after their detection 
( iPiro et al. 1998a|) . The ensuing discovery of GRB af- 
terglow emission in the X-ray, optical and radio bands 
has provided the distance for some GRBs and impor- 
tant information on the GRB remnants and their en - 



vironments, necessary to improve our understanding 
of the radiation mechanisms and the origin of GRBs. 
Time averaged spectra of GRBs seem to be consis- 
tent with the synchrotron shock model ( Tavani 1996 ), 
while Inverse Compton emission is expected to oper- 
ate at very early times ( Waxman 1997 ). 
Models of GRBs and afterglow emission are mainly 
based on mechanisms of dissipation of kinetic en- 
ergy of a relativistic expanding fireball (Meszaros and 



Rees 19971 , [Wijers et al. 1997| , |Vietri 1997D . There 
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are open issues concerning the production process of 
the prompt 7-ray radiation (e.g., internal or external 
shocks of the fireball material, Kobayashi et al. 1997), 
the time when the afterglow sets in as compared to 



the end of the prompt GRB emission (Sari 1997), the 
production process of the delayed radiation (e.g., ra- 
diative shock vs. adiabatic shock, see Sari, Piran and 
Narayan 1998 ). All these issues require observations. 

BeppoSAX offers the possibility to perform in a 
broad energy band (2-700 keV) temporal and spec- 
tral studies of the primary events from which after- 
glow emission has been detected. Indeed, when an 
event is simultaneously detected by the Gamma-Ray 
Burst Monitor (GRBM, 40-700 keV) ( [Frontera et~al 



1997|] |Feroci et al. 1997[ ) and by one of the two Wide 
Field Cameras (WFCs, 2-26 keV, Jager et al. 1997 ), 
it is possible to obtain, besides a precise localization 
(few arcmin radius error boxes), also X- and 7-ray 
spectra and time profiles of the GRBs. Therefore it 
is possible to study the spectral evolution of GRBs 
and relate the spectral and temporal properties of 
the main events with those of the associated X-ray 
afterglows. 

Results on the spectral evolution of GRBs in a broad 
band were reported by Strohmayer et al. 1998 , using 
the GRB detections obtained with the Gamma-Ray 
Burst Detector on board the Ginga satellite. How- 
ever, due to the lack of knowledge of the GRB di- 
rection and thus a poor knowledge of the instrument 
response function, those results were unavoidably of 
limited accuracy. 

We have already reported on the time averaged spec- 
tra of some GRBs observed with the BeppoSAX 
GRBM and WFCs ([Frontera et al. 1998bD . In this 
paper we will concentrate mainly on the evolution 
of the spectral properties of a sample of GRBs oc- 
curred between July 20, 1996 and April 25, 1998. 
For completeness we include in our sample also 
GRB960720 and GRB970228, the spectral evolution 
of which was already studied with results consistent 
with those obtained with the present analysis ( Pirc| 
et al. | 19984 [Frontera et al. 1998a| ). Results of com- 
parative spectral analysis of the GRBs detected after 
April 25, 1998, will be the subject of a future paper. 

2. INSTRUMENTATION AND GRB SAMPLE 

The GRBM consists of the 4 anti-coincidence 
shields of the Phoswich Detection System instrument 



(Frontera et al. 1997, Feroci et al. 1997, Amati et al 



1997[ ). Each shield is a CsI(Na) scintillator slab 10 
mm thick with a geometric area of 1136 cm^. Each 
slab is open to the sky, except along some directions, 
where material from other BeppoSAX instruments 
heavily decreases the transparency of the slab 7-ray 
entrance window. The GRBM detector operates in 



the 40-700 keV energy band. Two of the four GRBM 
units are co-aligned with the WFCs: unit 1 with 
WFC No. 1 and unit 3 with WFC No. 2. The data 
available from GRBM for spectral analysis include 
two 1 s ratemeters (40-700 keV and >100 keV) and 
240 channel spectra in the 40-700 keV band inte- 
grated over 128 s. For the GRB spectral evolution 
we use the 1 s ratemeters, but we check their consis- 
tency with the GRB time averaged spectra obtained 
from the 240 channel data( [Aniati et al. 1999| ). The 
energy resolution of GRBM units 1 and 3 with energy 
has been discussed by Amati et al. (1997) ; e.g., at 
279 keV (^o^Hg line) it is 20%. The on-axis effective 
area of the GRBM units 1 and 3 is 420 cm^ in the 
40-700 keV band and is 500 cm^ at 300 keV. 
The WFC instrument consists of two coded aperture 
cameras, each with a field of view of 40° x 40° full 
width at zero response and an angular resolution of 

5 arcmin. WFCs have an energy resolution ~ 20% at 

6 keV, are operated in normal mode with 31 channels 
in 2-26 keV and 0.5 ms time resolution. The on-axis 
effective area of WFCs No. 1 and 2 averaged in the 
2-26 keV energy band is 118 cm^ ( [lager et al. 1997 ). 
The follow-up observations of the GRB error boxes 
provided by the WFCs are performed with the Nar- 
row Field Instruments (NFIs) aboard BeppoSAX, 
that are orthogonal to both the WFC and GRBM 
axes. They include two imaging instruments (LEGS, 
0.1-10 keV, Parmar et al. 1997; MECS, 2-10 keV, 
Boella et al. 1997b ) and two direct-viewing detec- 
tors (HPGSPC, 3-60 keV, Manzo et al. 1997 ; PDS, 
15-300 keV, Frontera et al. 1997 ). 

The list of the BeppoSAX GRB events included in 
our sample is given in Table 1 along with some basic 
information: position error radius, offset angle of the 
GRBs with respect to the instruments (GRBM and 
WFC) axis, peak fluxes and time durations Tx and 
T^ in the X-ray (2-26 keV) and 7-ray (40-700 keV) 
band, respectively, time delay of the first NFI obser- 
vation from the primary event. Positive or uncer- 
tain (indicated by a question mark) X-ray, optical or 
radio afterglow detection is reported in the last col- 
umn. The GRB time durations are estimated from 
the 1 s ratemeters and give the time interval during 
which the GRB count rate is higher than the esti- 
mated background level by at least 2a. As can be 
seen, Tx is generally greater or at most equal to T-y. 

3. SPECTRAL ANALYSIS 

The light curves from the WFCs and GRBM in- 
struments are shown in fig. 1. Each light curve was 
divided into a given number of temporal sections (see 
fig. 1), and a spectral analysis in the 2-700 keV en- 
ergy band was performed for the average spectrum of 
each section. The duration of the sections was cho- 
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sen to be shorter during the rise of the burst in order 
to study the evolution of the primary event at the 
earhest times. The spectral analysis was performed 
f ollowing the same pro cedure used for GRB970228 
( [Frontera et al. 1998a ). The background level was 
subtracted from the GRB count rates (see fig. 1) as 
follows. For the GRBM spectra, this level was esti- 
mated using the count rates immediately before and 
after the GRBs. If the background is variable during 
the GRB, it is estimated by interpolation, using a 
quadratic function that fits 150 s count rate data 
before the burst and 150 s data after its end. For the 
WFC spectra and light curves, the background level 
is estimated using an equivalent section of the detec- 
tor area not illuminated by the burst or other known 
X-ray sources. We also checked the consistency of 
this background level with that obtained by using 
the data before and after the b urst. The respon se 
function of GRBM units 1 and 3 ( [Amati et al. 1997| ) is 
derived from on-ground calibrations and checked with 
the Crab Nebula that was cle arly detected using th e 
Earth occultation technique ( Guidorzi et al. 1998| ). 
The response function is currently known with an 
uncertainty of about 10% for incident photons with 
low offset angles with respect to the instrument axes, 
as it is the case of the GRBs in our sample (see Ta- 
ble 1). We added in quadrature this uncertainty to 
the Poissonian variance. 



We used XSPEC software package, issue 10 ( Arnaud 
1996|)| to deconvolve count rate spectra, assuming a 
given theoretical model as input function. In the 
following, all quoted errors correspond to 90% confi- 
dence level for each interesting parameter. 
A simultaneous fit to the time averaged WFC and 
GRBM spectra was performed by u sing as input 
mod el a photo-el ectrically absorbed ( [Morrison anc 
McCi immon 1983| ) smoothly broken power law (Banc 
et al. TOT 

given by: 



N{E) = Aexp{~aNH) 
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where a = cr{E) is the photo-electric cross-section 



of a gas with cosmic abundance ( Morrison and Mc- 
Cammon 1983), N// is the equivalent hydrogen col- 



umn density to the GRB, Tx and F^ are the power 
law low energy (below Eq) and high energy (above 
Eq) photon indices, respectively, and A is the nor- 
malization parameter. 



The derived best fit parameters, N//, Fx, F^, peak 
energy Ep of the logarithmic power (i^F(i^)) per pho- 
ton energy decade, are shown in Table 2 for each of 
the temporal sections. The value of Ep is given by 
Ep = £'o(2 + Fx), under the condition that F^ < —2. 
The reduced values obtained from the best fits 
were always acceptable (less than 1.1). We also de- 
rived lower limits on Ep when the i'F{i/) spectrum 
showed a bending at the 7-ray energies, but with 



F^ > 



-2. In this case, Ep was derived under the 



condition that the fit to the data was still acceptable, 
assuming a F^ slightly lower than —2 (= —2.1). The 
lower limits on Ep reported in Table 2 were thus ob- 
tained from the lower limits on Eq corresponding to 
a 90% confidence level. 

When the N// values were not constrained from the 
data, they were fixed at the Galactic values along the 
GRB directions, except for GRB980329, for which we 
adopted the value derived from the late afterglow ob- 



servation (in 't Zand et al. 1998). 
Also reported in Table 2 are the spectral and tem- 
poral properties of the associated X-ray afterglows 
in the 0.1-10 keV energy band as derived from the 
TOO observations with the NFIs. The two spec- 
tral parameters (photon index F^ and column den- 
sity N//) were derived assuming a photo-electrically 
absorbed ( [Morrison and McCammon 1983 ) power law 
model with photon index F^, while the parameter 5 
is the index of the fading law I(t)oc t~^, which best 
fits the 2-10 keV light curve. 

4. RESULTS 

For three GRBs (GRB960720, GRB970508, 
GRB980425) the 7-ray emission starts earlier than 
the X-ray emission, while in the other GRBs 7-rays 
and X-rays rise simultaneously (fig. 1). Duration 
and shape of the GRB time profiles change from 
one GRB to the other, even if some similarity is ob- 
served between GRB980329 and GRB980425, while 
GRB960720 and GRB970508 exhibit strikingly sim- 
ilar light curves. No X-ray precursor activity is de- 
tected, similar to that seen in other GRBs by the 
Gamma Burst Detector (GBD) on-board Gin ga ([Mu- 
rakami et al. 199l| ) or WATCH/Granat ( [Castro- 



Tirado et al. 19941) . 



Time duration and shape of the GRB time pro- 
files change from one GRB to the other, even if 
some similarity is observed between GRB980329 and 
GRB980425. 

From the low values (<1-1) of the reduced can 
infer that the GRB spectra of our sample can be de- 
scribed by a Band law (see eq. 1), even if, in some 
temporal sections (in the case of GRB970402 in both 
sections), they are single power laws. GRB spectra 
evolve with energy and with time (see Table 2). The 
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spectra are generally harder at low energies or they 
exhibit the same slope at high and low energies. We 
do not find bursts with low energy excesses with re- 
spect to the Band law (see eq. 1) as claimed by other 
authors ( Preece et al. 199^ , Strohmayer et al. 199§| ). 
In the first few seconds, during the GRB rise (fig. 1) 
we find, for 4 GRBs, that the low-energy photon in- 
dex Tx is significantly larger than -2/3, limit slope 
of the low energy tail of the instantaneous energy 
spectrum of an optically thin synchrotron emission 
(see Section ^]^). This feature was already noticed 
for GRB960720 (|Piro et al. 1998a| ). In one case 
(GRB970111) it holds for the entire GRB duration. 
The spectra are generally softer during the tail of the 
GRBs. A remarkable exception is observed in the 
spectrum of GRB970228, that shows a softening dur- 
ing the first pulse, but it becomes harder during the 
subsequent minor pulses of the event tail (see fig. 1). 
This feature was already discussed by Frontera et al. 
(1998a) . 

If we compare the GRB photon index Tx during the 
GRB tail with the photon index Ta of the power 
law spectrum of the associated X-ray afterglow, when 
available, it can be seen that the afterglow spectrum 
is generally softer than the GRB tail spectrum. 
Generally the absorption column densities are con- 
sistent with the Galactic values along the GRB di- 
rections, except for GRB980329. In this case, the 
Nh values are higher than the Galactic value (0.94 x 
10^^ cm~-^) in the temporal sections from A to F, 
while in the last two temporal sections (G and H), 
we are not capable to constrain the N/f value, that 
was frozen to the value estimated during the late 
afterglow observation ( [in 't Zand et al. 1998| ) , 
but was also consistent with the Galactic value. A 
much higher hydrogen column density (about Ix 
10^'^ cm~^) is specially apparent during the section 
E corresponding to the early decay of the primary 
event (see fig. 1 and Table 2). 

The ly'F^u) spectrum in each of the temporal sections 
in which we subdivided the time profiles of the GRBs 
in our sample is shown in fig. 2. The evolution 
of the peak energy Ep is apparent for most of the 
GRBs studied. As an example, GRB970228 shows 
a peak energy quickly crossing the 2-700 keV pass- 
band during the first peak, while the peak energy of 
GRB970111 is still within that energy range until the 
end of the event. There is no evidence of softening 
in the spectra of GRB970402, but from the slope of 
the z/F(j/) spectrum (see fig. 2), it is clear that Ep 
is above 700 keV. GRB980425, presumably associ- 
ated with the supernova SN 1998bw ( [Galama et al. 
1998lj), does not show special features with respect 
to the other GRBs in our samples: its peak energy 
decreases with time from the GRB onset. 



The GRB970111 behavior deserves particular atten- 
tion. The X-ray spectral shape is unchanged, while 
the 7-ray spectrum steepens monotonically. This in- 
dicates a progressive decrease of Ep, which however 
almost stabilizes at 30-40 keV starting about 23 s af- 
ter the GRB onset. This result is better apparent in 
fig. 3, that shows the time behavior of Tx as a func- 
tion of for all GRBs considered. When the peak 
energy is above or below the 2-700 keV passband, we 
would expect to measure approximately similar pho- 



ton indices Tx and F^ (see Section 5^) and the value 
of (Fx, F^) should be along the diagonal (dashed line 
in fig. 3). On the contrary. Fx and F^ could be 
markedly different (with Fx close to the limit value - 
2/3 for an optically thin synchrotron spectrum) when 
Ep is in the 2-700 keV passband. In this case (Fx, 
F^) should lie out of the diagonal. If Ep, during the 
primary event, sweeps the above passband, the plot 
of (Fx, F^) should start on the diagonal, go away 
from it and go back on the diagonal. In the case of 
GRB970111 Ep enters the 2-700 keV passband soon 
after the onset, but does not completely cross it un- 
til the end of the GRB. In some cases (GRB960720, 
GRB970228) it sweeps 2-700 keV passband, while in 
other cases (GRB971214, GRB980329, GRB980425) 
the peak energy is still in the passband at the end of 
the primary event. In the case of GRB970402, from 
the rising of its I'Fi' spectrum (fig. 2), the peak en- 
ergy is clearly above 700 keV, while for GRB970508 
there is a hint that Ep crosses the 2-700 keV pass- 
band. 

In only one case (GRB971214) we observe a harden- 
ing in the afterglow spectrum with respect to that 
observed during the prompt emission as expected in 
the model discussed by Waxman (1997) . 

5. DISCUSSION 

We discuss our results in the light of current the- 
ories on GRBs in order to contribute to solve many 
open issues concerning the GRB phenomenon. 

5.1. Test of the Synchrotron shock model 

As discussed by several authors (see, e.g., Katz 
1994, Tavani 1996, Sari k Piran 1997b), the electro- 
magnetic radiation from GRBs and their afterglows is 
likely due to synchrotron radiation from accelerated 
particles in internal or external shocks of a relativis- 
tically expanding fireball (Synchrotron Shock Model, 
SSM). The observed hard-to-soft change of the spec- 
trum with energy and with time from the GRB onset 
is in agreement with this picture ( permer, Bottcher 
and Chiang 19991 ). 

Assuming that the emission region is optically thin, 
the expected power law index of the vT{v) spectrum 
below the peak energy Ep is expected to be in the 
range from 1/2 to 4/3 (corresponding to -3/2 to - 
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2/3, respectively, in the photon spectrum), where the 
first index is obtained when coohng of the particle 
distribution during the GRB is also taken into ac- 



. The other possibility is that intrinsic absorption 
takes place at the GRB site. If this is the case, in 
the context of the fireball shock model ( |Sari, Piran 



count, while the latter index is reached in the case and Narayan 1998), N// is given by E/ (7onipCr)di 



of an instantaneous spectrum or for monochromatic 
particles ( Cohen et al. 19^ ). Notice that the limit 
of -2/3 is independent of the electron distribution or 
uniformity of the associated magnetic field. 
As discussed by Tavani (1996) , the Band law is a 
good approximation of the spectrum expected by the 
thin SSM and our results confirm this fact. However, 
using this simple model, we find that, for 50% of the 
GRBs (see fig. 4), Fx is above the expected limit pho- 
ton index (-2/3) in the earliest time intervals of the 
light curve. In the case of the intense GRB970111, 
this holds almost for the entire GRB duration. Crider 
et al. (1997) found a similar violation in the 30- 
1800 keV spectra of GRBs observed by BATSE, while 
Strohmayer et al. (1998) and Preece et al. (1998) 
find this violation in time averaged spectra of GRBs 
detected with the GBD/Ginga and BATSE/CGRO 
experiments, respectively. 

Our data confirm that some additional mecha- 
nism is certainly active for some GRBs, that mod- 
ifies the optically thin synchrotron spectrum at 
very early times. This could be synchrotron self- 
absorption ([Papathanassiou 19991), Compton up - 
scattering of low-energy ph otons ([Liang et al. 1997| ) 
or plasma physics effects ( p?avani 1999| ). Assuming 
single Compton up-scattering by highly relativistic 
electrons with random direction, the limit photon in- 
dex of the low-energy spectrum is zero ( Rybicki and| 
Light man 1979). We find that the low-energy pho- 
ton index of all GRBs in our sample, except that of 
GRB970111, is consistent with this limit value. 
The approximate stabilization of Ep at ~35 keV for 
GRB970111, corresponding to a decrease of the power 
law photon index F^ will be discussed in Section 5.3.3 . 
For the other GRBs, the change of the photon index 
observed starting from the GRB tail up to the late 
X-ray afterglow observations not as much relevant. 
The high Njj observed during the prompt emission 
of GRB980329 and its decrease during the tail of 
the event can be interpreted in two ways. It might 
be a consequence of variable absorption of external 
intervening material due to its progressive photo- 
ioniziation as the GRB evolves ( [Bottcher et al. 1998| ) . 
In this case Bottcher et al. (1998) also expect fluores- 
cence features and/or K edges that are not appar- 
ent in the prompt and delayed emission ( in 't Zand| 



et al. 199'^ ). However these features could be too 



weak to be detectable with our instruments and/or 
could be not within the WFC or NFI passband, due 
to the redshift. The redshift of GRB980329 is still 
unknown (see a possible estimate in Fruchter 1999) 



vided by Anr'^, where E is the energy of the shock, 
7o is initial Lorentz factor of the shocked material, 
uip is the proton mass, c is the light velocity ad r is 
the distance from the explosion center to the shock 
dissipation radius. With typical values for the en- 
ergy E (~ 10^^ erg) and 70 (~ 100), the measured 
Nh gives a value of r~ 10^^ cm, which is inconsistent 
with external shocks that are expected only at values 

of R-^ 10^^ 10^^ cm, when the N// would be lower 

by a factor of 10^ — 10^. Thus our results points rather 
to internal shocks accompanied by rapid expansion of 
a fireball. 

5.2. Earliest X-ray afterglow emission 

Table 3 shows the X-ray (2-10 keV) and 7-ray (40- 
700 keV) fluences, Sx and S^, of the GRBs in our 
sample as derived from our spectral analysis. For 
comparison, when available, also the 2-10 keV flu- 
ence Sa of the afterglow emission is given. The latter 
was estimated from the fading laws (see Table 2) in- 
tegrated over the time interval from the GRB end up 
to 10^ s, when the X-ray afterglow emission level is 
generally already negligible for a power slope 5 of the 
fading law greater than 1, as in our case. The ratio 
Sx/S-y changes from event to event and ranges from 
1% to about 40%. It does not seem that detection of 
X-ray afterglow emission is related to the amplitude 
of the above ratio. Also the ratio Sa/S^ changes from 
a GRB to another, ranges from 0.4% to about 20% 
and does not appear to be strictly related with the 
Sx/S-y value, as also evident from the Sa/Sx ratio. 
This ratio also shows that the fluence of the X-ray 
afterglow emission is, within a factor 2, of the same 
order of that measured, in the same energy band, 
during the primary event. 

In the scenario of the fireball model, in order to 
explain the complex time profiles of GRBs, it was 
shown (Kobayashi et al. 1997, Sari and Piran 1997b| ) 
that the GRB main event could be due to inter- 
nal shocks, although for simple fast-rise/exponential- 
decay (FRED) events also external shocks can be ac- 
cepted (Permer, Bottcher and Chiang 1999). In- 



stead, there is general consensus that the GRB late 
afterglow emission is likely due to the external shocks 



propagatin g in the interstellar medium ( Meszaros and 



Rees 1997 , 3ari, Piran and Narayan 1998| ). In this 
scenario, the two phenomena can evolve in different 
ways and can involve energetics not necessarily corre- 
lated. However it has been suggested ( pari 1997 , Sari 



and Piran 1999b| ) that the early afterglow could start 



few dozens of seconds after the burst. We tested this 
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prediction for the GRBs in our sample. We have al- 
ready shown (see, e.g., Costa et al. 1998 and Piro 
et al. 1998 ) that the extrapolation of the 2-10 keV 
afterglow fading law back to the time of the GRB is 
in agreement with the 2-10 keV flux measured during 
the last portion of the event. This should suggest that 
at least the tail of the GRBs could be due to after- 
glow emission, while at the beginning of the primary 
event, the contribution to the GRB emission from the 
afterglow is negligible. We adopt the following proce- 
dure to perform the test. It is reasonable to assume 
that the intensity of the X-ray afterglow depends on 
the spectral index Ta and temporal index 6 according 
to the following law: 



I{E,t) oc E^H 



(1) 



With this assumption, the more general expression 
of the ratio between the 2-10 keV fluence in the time 
intervals (ti, 12) and (t2, 00), with t2 >ti, is given by 

R{h,t2) = Sih,t2)/S{t2,00) = i^(fi,f2)[(tl/t2)"'^+'- 

(2) 

where S(ti,t2) is the 2-10 keV fluence in the (ti,t2) 
interval, S(t2,cxD) is the fluence in the interval from t2 
to a time t ^ t2 where the afterglow intensity has de- 
creased to a negligible value, Fi and T2 are the aver- 
age photon indices of the afterglow power law spectral 
emission in (ti, t2) and (t2, 00), respectively. Notice 
that parameter i^(ri,r2) must be greater than 1 if 
the spectrum becomes softer with time as observed 
(see Table 2). 

If we express ti and t2 in terms of the time duration 
T of each GRB, that is, ti = fiT and t2 =f2T, with 
fi <f2, the previous equation becomes 



Rifl,f2) = K(fi,f2)[(/l//2) 



-5+1 



1] (3) 



The ratio R(fi,f2) for each GRB of our sample was 
evaluated assuming f2 = 1 and different values of fi . 
The result obtained for fi = 0.01 and fi = 0.63 is 
shown in fig. 5. The above law gives a good fit to 
the data {xl = 1-04, dof=4) in the case fi = 0.63, 
but an unacceptable fit (x^ = 2.5, 4 dof) in the case 
fi = 0.01 when the total GRB fluence is used. We 
also checked the behavior of the above ratio in the 
case fi = 0.40, obtaining a reduced Xu — In 
fig. 5 (bottom panel) we show the best fit curve, with 
K = 1.66ib0.36, in the former case. For comparison, 
the curve with K=l, under the assumption of no spec- 
tral evolution from (ti, t2) to (t2, 00), is also shown. 
In fig. 5 (top panel), we report the expected curve 
with the minimum value of K allowed (K=l) in the 
case fi = 0.01. There is no agreement between data 
and expected dependence of R(0.01,l) on 5. 
This result shows that roughly the second half of 



the GRB X-ray light curve is consistent with after- 
glow emission. Likely GRB time profiles are the su- 
perposition of two components, one due to the pro- 
cess that produces the main event and the other due 
to the delayed emission. A possible scenario envis- 
ages the former component as due to internal shocks 
and the latter due an external shock (see, e.g., Sari 
(1997)). If the emission due to internal shocks has a 
short duration, the afterglow can be separated from 
the internal component, that gives rise to the ear- 
liest emission of a GRB. This could be the case of 



GRB970228 (|Costa et al. 1997| , [Frontera et al. 1998aD , 
in which the first peak (see fig. 1) could be due to an 
internal shock, while the late emission could be due 
to external shocks and thus coincide with the after- 
glow emission. 

According to the internal-external shock model for 
the GRB, from the above results, it is possible to 
evaluate the initial Lorentz factor 70 of the shocked 
fluid ( pari and Piran 1999a| , |Sari and Piran 1999b| ): 



-1] 



70 



240S 



1/8 1 



52 



n^ 



\t/W s) 



-3/8 



(4) 



where t is the onset time of the afterglow emission, 
E52 is the afterglow energy released in the shock in 
units of 10^^ erg and ni is the number density of the 
ambient medium in units of cm~^. If the early after- 
glow starts at about 50% of the GRB duration, for 
those GRBs for which it is possible to make an es- 
timate of E52 (see Table 3), we obtain the values of 
7o reported in Table 4, assuming a constant ambient 
medium density ni = 1. The luminosity distance 
was evaluated from the redshift values reported in 
Table 3, assuming a standard Friedman cosmology 
with Hq = 65 km/s/Mpc and 17 = 0.3. The energy 
released was assumed to be isotropic. The value of 70 
is of the order of 150 for all GRBs in our sample with 
known distance, except for GRB980425. For this pe- 
culiar GRB, we have assumed its coincidence with the 



type Ic supernova SN1998bw ( [Galama et al. 1998b|) . 
In this case. Plan et al. 1999b estimate an afterglow 
emission with fading law index 6 > 1.3. This upper 
limit was obtained assuming that the afterglow starts 
during the GRB tail as discussed above, and consid- 
ering only the first TOO measurement as main contri- 
bution to the afterglow radiation, while the radiation 
detected about 6 months after the main event was 
assumed to be mainly due to thermal emission from 
the supernova. The much lower value of 70 (< 50) 
obtained for this GRB shows that the association of 
GRB980425 with SN1998bw implies a deviation from 
the extreme relativistic conditions expected for GRBs 
in the context of the fireball scenario and thus that it 
could be a member of another class of GRBs. Assum- 
ing that this event is not associated with SN1998bw, 
but with a fireball with redshift z ~ 1, we could ob- 
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tain a value of 70 consistent with the values found for 
the other GRBs. 

5.3. Hydrodynamical evolution of a shock 

An observational test of the hydrodynamical evo- 
lution of GRB primary events and afterglows is of 
key importance to understand the GRB phenomenon. 
The afterglows evolution resulting from forward 
shocks with the interstellar medium (ISM) was the- 



oretically investigated by several authors (f 


Sari, Piran 


and Narayan 1998, 


Meszaros and Rees 1997 


, pVIeszaros 



and Rees 1999 ), while the GRB evolution is less clear, 
except the case of FRED (fast rise, exponential de- 
cay) events ( permer, Bottcher and Chiang 1999| ). 
Here we test the expectations of the model by Sari, 
Piran and Narayan (1998) making use of our spec- 
tral results. This model considers a spherical shock 
that propagates in a surrounding medium (interstel- 
lar medium, ISM) of constant density. The electrons 
accelerated in the shock are assumed to have a power 
law distribution of Lorentz factor 7e (N(7e) oc 7"^) 
where 7e = Ee/irieC^ > 7m and electron distribution 
index p > 2. Two extreme cases of hydrodynamical 
evolution of the shock are considered: a fully adia- 
batic and a fully radiative evolution. In the case of an 
adiabatic evolution the internal energy of the shock 
is constant, while in the other case it varies with the 
bulk Lorentz factor 7. The time dependence of 7 is 
different in the two regimes. 

5.3.1. Late GRB spectrum vs. afterglow decay 

If the GRB tail is due to afterglow emission and 
thus to an external shock, we should also expect a 
correlation between spectral properties of the GRB 
tail and power-law index 6 of the fading afterglow 
emission. The model by Sari, Piran and Narayan 
(1998) predicts the occurrence of two breaks in the 
synchrotron spectrum, one at energy E^ correspond- 
ing to the lowest energy of the injected electrons 
with Lorentz factor 7^, and the other at energy Ec 
(cooling break), corresponding to the electron energy 
with Lorentz factor 7c, above which cooling for syn- 
chrotron radiation is significant. Figure 6 shows the 
expected evolution of the i^F^v) spectrum for fully 
adiabatic and fully radiative cooling of the shocks, 
for typical values of the model parameters (see val- 
ues in the figure caption). Due to the different vari- 
ation with time of Ec and Em (see Sect . |5. 3.2 ), at 
early times, Em >Ec (fast cooling), while at late 
times, Em <Ec (slow cooling). The passage from 
fast cooling to slow cooling occurs at a time to when 
Em = Ec = Eq. It was shown ( |Sari, Piran and| 
Narayan 1998| ) (see also fig. 6) that the expected 
power law photon index of the afterglow spectrum 
is given by -3/2 for Ec <E<Em and by -p/2-1 in the 



high energy tail (E higher than Em and Ec), inde- 
pendently of the cooling mechanism of the shock (ra- 
diative or adiabatic) . From Fig. 2 and its comparison 
with fig. 6, we can see that, in the case of GRB970228 
and GRB970402, the high energy photon index Vi^ of 
the last section of the time profiles is consistent with 
-3/2, while for the other GRBs in our sample it is 
likely related with the high energy tail of the asso- 
ciated afterglow. We make this assumption setting, 
for all GRBs in our sample except GRB970228 and 
GRB970402, F^ = -p/2-1. 

The fading law of the 2-10 keV afterglow emission 
(I(t)(xt~'^) provides a further relation between p and 
S. Indeed the model by Sari, Piran and Narayan 
(1998) expects that the index 6 is related to p through 
6 = 3p/4 — l/2in the case of an adiabatic cooling 
of the shock when both Em and Ec are below the 
2-700 keV passband at the epoch t of the afterglow 
observation. This is likely our case. The same model 
expects that, in the case of a radiative cooling of the 
shock (in this case the shock becomes adiabatic at 
time to when Em = Ec), 6 = 6p/7 — 2/7. Com- 
bining all the above relations, we found the following 
dependence of F^^ on 6: Vi^ = —4/3 — 2/36 for an 
adiabatic shock, and F^ = —7/6 — 7/125, for a ra- 
diative shock. 

The measured values of Fi^ as a function of 6 are 
shown in fig. 7a. Superposed to the data are the 
curves derived from the above relations. For compari- 
son, we show in Fig. 7b, for all GRBs in our sample for 
which the afterglow spectrum is known, the power- 
law photon index F^ of the X-ray afterglow (see Ta- 
ble 2) function of the corresponding index 6 of 
the fading law in the same energy band. Also shown 
in the figure are the expected F^ vs. 5 relationships 
pari, Piran and Narayan 1998| ): F^ = -7/6-7/12(5 
for radiative cooling (t<to), while for an adiabatic 
cooling, F„ = -2/35 - 1 or F^ = -2/36 - 4/3, de- 
pending on the observing time t after the primary 
event, if tm <t<tc or t>tc, respectively. The param- 
eters tm and tc are the times at which Em and Ec 
sweep the energy passband of the MFCS detectors, 
respectively (see fig. 6). 

From fig. 7a, it appears that the data are consis- 
tent with the expected dependence of F^ on 5, but 
their uncertainties do not allow us to discriminate 
among different cooling models. The peculiarities of 
GRB970111 and GRB980425 are also apparent from 
this plot. Assuming that the afterglows associated 
with these GRBs satisfy one of the two plotted rela- 
tionships, the X-ray afterglow of GRB970111 should 
have exhibited a fading law with power-law index 
lower than -2, that is consistent with its non detection 
with BeppoSAX 17hrs after the primary event. The 
same condition applied to the GRB980425 afterglow 
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implies a much lower 6 than -2. 

Also the measured Fq vs. 5 are consistent with the 
expected relationships between these two quantities 
derived from the synchrotron shock model( pari, Pi- 
ran and Narayan 1998|), but, also in this case, due 



to the uncertainties in these parameters, we cannot 
discriminate among different cooling laws. 

The above discussed 5 vs. p and T vs. p relation- 
ships have also permitted the estimate of the index of 
the electron distribution for 5 GRBs in our sample. 
Given the lower uncertainty in the value of 6 (see Ta- 
ble 2), we report in Table 4 the values of p derived 
from 5. We assume an adiabatic cooling (see sect. 
5. 3. 2D in the two cases above considered (tm <t<tc, 
and t>tc). The corresponding values of p derived 
from Ta, due to the larger uncertainty in this param- 
eter (see Table 2) are consistent with the p values 
reported in Table 4. These change from one GRB to 
the other and rage from 2.1 to 3.1. It is remarkable 
that for the two GRBs (971214 and 980329) for which 
it was possible to derive the values of p from the spec- 
tral index Vi^ of the prompt emission, we found val- 
ues of p also consistent with those reported in Table 4 
(3.2±1 for GRB971214 and 2.6±0.6 for GRB980329). 

5.3.2. Temporal evolution of the peak energy Ep 

As shown in fig. 6, the model by Sari, Piran and 
Narayan (1998) predicts a temporal evolution of the 
i'F{i/) spectrum of the GRB afterglow. A spectral 
evolution is also expected during the primary event, 
even if the time behavior can be more complex: it can 
change from one GRB to the other and can depend 
on the assumed model (e.g., external shocks, Dermer, 
Bottcher and Chiang 1999 vs. internal shocks. Sari 
1997 ). We observed both spectral evolutions (see Ta- 
ble 2 and fig. 2): on the basis of the results discussed 
in the previous sections, the spectra obtained in the 
first half of the GRB duration are mainly due to the 
primary event, while those obtained later are likely 



due to afterglow emission (see Sect. 5. 3.1 ). 
To characterize the spectral evolution, we consider 
the time behavior of the peak energy Ep of the z^F(i/) 
spectrum. The model by Sari, Piran and Narayan 
(1998) predicts the occurrence of a maximum in the 
z^F(i/) spectrum in correspondence of one of the two 
breaks in the synchrotron spectrum, at energies E^ 
and Ec, respectively (see Sect.^^). These parameters 
decrease with time in a different manner and accord- 
ing to the assumed cooling type of the shock (adia- 
batic or radiative) . At early times from the afterglow 
onset (t<to) the peak energy Ep is coincident with 
Em that decreases as t^^/^ for adiabatic shock and 
as t""'^^/'^ for a radiative shock, while at late times 
(t>to) it is coincident with Ec that decreases as t~^/'^ 
for an adiabatic shock and as t"^/'' for a radiative 



shock (see fig. 6). As discussed above, at the time 
to Em coincides with Ec (and with Ep). This time is 
different for an adiabatic or a radiative cooling (see, 
e.g., fig. 6). Notice that these temporal dependences 
of Em and Ec hold during the decay phase of the 
afterglow, while at the beginning of the fireball inter- 
action with the ISM one expects a shallower decay of 
the break energies. 

Thus the observed time behavior of Ep can constrain 
the cooling mechanism of the shock and provide the 
value of important physical parameters, like the frac- 
tion Ec of shock energy that goes into the electrons 
and the fraction e b that goes in magnetic energy den- 
sity ( Sari, Piran and Narayan 1998| ). 
Figure 8 shows the behaviour of Ep with respect to 
time for all GRBs in our sample. The values are 
not corrected for GRB redshift, even when this is 
known. We have verified that this correction is not 
critical for our conclusions. For the GRBs with de- 
tected late afterglow, we also show in fig. 8, when 
available, the values of Ep obtained from the litera- 
ture (GRB970228, Masetti et al. 1999 ; GRB970508, 
Galama et al. 1998a ; GRB980329, Palazzi et al. 1998 
; GRB971214, Dal Fiume et al. 1999 , Ramaprakash 
et al. 1998 ). Superposed to the data is also the ex- 
pected time behavior of E^ for short times from GRB 
onset in both the radiative and adiabatic conditions 
of the shock. In addition we show the behavior of 
Ec when afterglow emission is detected and the spec- 
trum is known. The E^ curves are normalized to the 
earliest measured values of Ep, while the Ec curve 
is normalized at the value of Ep obtained from the 
late afterglow spectrum. In fig. 8 the time at which 
the afterglow is expected to have already started (see 
Sect j5.2| ) is also shown (vertical dashed line). 

For most of the GRBs in our sample, only lower 
limits of Ep are available in the earliest part of the 
event. In spite of that, it is apparent that Ep de- 
creases with time more quickly than expected in the 
case of a fast cooling (adiabatic or radiative) of an 
external forward shock. This is another evidence for 
internal shocks during the first part of GRBs. 
Another important feature of the data is that, for 
GRB970228 and GRB970508, the extrapolation of 
the fast cooling curves at late times is clearly incon- 
sistent with the value of Ep measured during the af- 
terglow about 10 days after the primary event. This 
means that a transition from a fast cooling to a slow 
cooling phase of the shock has occurred in the mean- 
time. For GRB970508 we report in fig. 8, in addition 
to the late value of Ep =Ec, the estimated value of E^ 



( palama et al. 1998a ). This break fits nicely with 
the adiabatic evolution curve, confirming the passage 
of the shock from fast to slow cooling. The backward 
extrapolation of the expected behavior of Ep normal- 
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ized to the last data point allows us to derive the 
transition time to from fast cooling to slow cooling. 
The to values obtained are reported in Table 4, along 
with the rough estimate of the same parameter for 
other two GRBs (GRB971214 and GRB980309), for 
which the transition time to is more uncertain. 
The Ep behavior alone does not permit us to clearly 
establish whether at early times the shock cooling is 
adiabatic or radiative. We obtained this information 
by comparing the maximum of the energy spectra 
Fmarr at early time (during the latest section of GRB 
time profiles) with that at late time (during late af- 
terglow observations). The results are shown in Ta- 
ble 4. In the case of GRB970508 and GRB980329, 
Fmax (early) and FmaxO-SLte) are comparable, while in 
the case of GRB971214, F^aa; (early) is a factor from 
about 4 to 18 higher than Fmax(late). This range is 
due to the uncertainty in the true value of Fmax (late) . 
The two values reported in Table 4 are those given 
by Ramaprakash et al. (1998) and Dal Fiume et al. 
(1999) , respectively. The difference is due to the dif- 
ferent absorption by dust assumed by these authors 
to correct the measured magnitudes. Dal Fiume et 
al. (1999) assumes the extinction typical of starburst 
galaxies and take into account the redshift of the mea- 
sured radiation, while Ramaprakash et al. (1998) as- 
sume an exponential absorption with optical depth 
from extinction oc 1/A, where A is the wavelength of 
the observed radiation. 

As a result, at early times, for GRB970508 and 



GRBB80329, we see an adiabatic evolution of the 



shock, while for GRB971214, the shock cooling is 
close to be adiabatic if we assume the highest value of 
Fmax (late). As discussed by Dal Fiume et al. (1999) 
, there are reasons to prefer their correction, which 
implies that the shock is almost adiabatic. 
Also the values of the shock parameters eg and €b are 
reported in Table 4. In the case of GRB970228, we 
are not able to estimate Fmax (late) from the avail- 



able data (iMasetti et al. 1999|) . For this GRB we 
report the values of these parameters for both a ra- 
diative and an adiabatic fast cooling. It is noter- 
worthy that the derived eB(late) and (early) are 
roughly in agreement and that they are lower than Cg- 
This is expected when the contribution from Comp- 
ton up-scattering cannot be neglected. The value of 
to derived by us for GRB970508 is consistent with the 
value (about 700 s) derived by Galama et al. (1998a), 
who also noticed that an adiabatic model fits better 
the Idte afterglow data. 



5.3.3. The peculiar behavior of GRB970111 

We already have discussed in Sections 4 and 5.1 
some peculiarities of GRB970111. This GRB stands 
out for its low energy spectral slope, that is not con- 



sistent with a thin synchrotron model and with In- 
verse Compton, making self-absorption a possible ex- 
planation. GRB970111 is also peculiar for its Ep time 
behavior: a transition in the rate of Ep decrease is 
observed in correspondence of the expected starting 
time of the afterglow emission (see fig. 8). However 
this break cannot be interpreted as passage of the 
cooling break Ec through the 2-700 keV passband, 
since the change of the spectral slope above and below 
Ec is predicted to be 1/2 and we see a much greater 
change. It therefore must be E^ which should how- 
ever decrease as t~'^/^ in the adiabatic case or faster 
(as t~^^/^) in the radiative case. Instead the slope of 
Ep is consistent with t~^/^ that is not in agreement 
with the decay phase of the afterglow. 
The evolution of the shock does not appear to be 
adiabatic. The non adiabatic character is apparent 
from the time behavior of Fmax ■ In the second half of 
the GRB time profile, the energy spectrum shows an 
intensity maximum that decreases with time accord- 
ing to a power law with index 2.6, while a constant 
intensity is expected for an adiabatic cooling. The 
measured slope is also inconsistent with a radiative 
cooling, in which case a much lower slope (-3/7) is 
expected for the power law time decay of F^ax^ 

The 

extrapolation of the time behavior of Ep and F^ax 
at the time of the first BeppoSAX observation yields 
X-ray fluxes that are well below the upper limit given 
by Feroci et al. (1998) . Also the upper limits of the 



afterglow given in the optical band (Gorosabel et al 
1998| ) are well above this extrapolation. 



A possible intepretation of this strange Ep behavior 
is that we are observing two successive electron accel- 
eration episodes, giving rise to the first and the sec- 
ond pulses (see 7-ray time profile in fig. 1). During 
the first pulse, Ep is supposedly decaying but never 
enters the 2-700 keV passband (as soon as injection 
stops the peak can sweep this energy band faster than 
allowed to observe with our time integration of the 
spectra), while during the second pulse the electron 
distribution has a long cooling timescale that makes 
the peak almost stable at around 35 keV. This second 
episode (that could also be originated by an external 
shock) should show a roughly constant bulk Lorenz 
factor that later decays and gives all the power de- 
pendences on the lightcurve. Following this interpre- 
tation, the electron power law dependence on F^ is 
p= 1 -2r^ for the early portions of the second pulse 
and p= -2(r^+ l) for the later ones ( |Sari, Piran and 
Narayan 1998| ). From the measured values of F^ we 
find, for the second pulse of the GRB, p~4. This 
electron index can account for the non detection of 
X-ray/optical afterglow emission about 13 hrs after 
the primary event. 
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6. CONCLUSIONS 

Several conclusions can be drawn from the above 
results. The optically thin synchrotron model (e.g. 
Tavani 1996 ) appears to be consistent with the GRB 
spectra reported in this paper. It can describe most 
of them and their evolution with time and energy, 
but not the low energy spectra at earliest times after 
the GRB onset, where there is a deviation from its 
expectations. Some additional ingredient, probably 
self-absorption or Compton up-scattering, is required 
to explain the low energy spectra at early times. In 
the case of GRB970111, the deviation is the largest 
and lasts for the entire GRB duration. 
In one case, GRB980329, a relevant and variable ab- 
sorption has been detected during the prompt emis- 
sion. This feature appears to be consistent with vari- 
able self-absorption of the X-ray emission produced 
in internal shocks of a relativistic expanding fireball, 
but absorption from external material that progres- 
sively becomes photo-ionized by the GRB radiation 
cannot be excluded. 

The fluence of the second half of GRB time profiles 
is consistent with that expected from afterglow emis- 
sion. With this assumption, in the context of the 
fireball shock model (e.g. Sari, Piran and Narayan 
1998 ), we have derived the initial Lorentz factor 
7o of the shocked material for 5 GRBs in our sam- 
ple, for which the redshift has been determined. The 
value of 7o is about 150 for all these events, except 
for GRB980425 if this GRB is associated with the 
supernova SN1998bw. In this case we find an upper 
limit 7o < 50, that is lower than the typical values as- 
sumed in relativistic expanding fireballs. This GRB 
could be a member of another class of events. 
We have tested the presence of a correlation between 
photon index of the late energy spectra of the GRB 
prompt emission and index of the fading law of the 
associated afterglow. We find that this correlation ex- 
ists and it is consistent with the hydrodynamical evo- 
lution of an external shock as discussed in the model 
by Sari, Piran and Narayan (1998). However the un- 
certainty in the data does not allow the discrimina- 
tion between different types of evolution (adiabatic 
vs. radiative, fast cooling vs. slow cooling). 
The shock cooling properties have been investigated 
from the temporal behavior of the peak energy of the 
i^F(z/) spectrum. We have still assumed the external 



shock model by Sari, Piran and Narayan (1998) . We 
find that the values of Ep measured in the first half 
of the GRB time profile have a time behavior which 
is not consistent with an external shock and confirm 
that internal shocks have a key role in determining 
the GRB spectral evolution at earliest times. 
Comparing the values of in the second half of the 
event with those extimated at late times (~ 10^ s 
from the GRB onset) when available, we find that for 
most of these events, a transition from a fast cooling 
to a slow cooling has taken place. Adiabatic evolu- 
tion seems to be preferred. In the case of GRB971214, 
the uncertainty in the peak fiux of the afterglow spec- 
trum does not allow the evolution to be determined, 
although an almost adiabatic evolution is preferred. 
We also have determined the fractions of the shock 
energy that go in electrons and magnetic energy, and 
the index p of the electron distribution accelerated 
in the shock. We find that p changes from one GRB 
to the other and ranges from about 2.1 to 3.1, for an 
adiabatic shock cooling. 

The case of GRB970111 merits particular considera- 
tion. Prom this GRB no afterglow emission was de- 
tected about 17 hrs after the primary event. The 
GRB time profile in the 40-700 keV energy band 
shows two main pulses, with the second one that oc- 
curs in the second half of the GRB duration. We 
find that the time behavior of the peak energy of the 
z/F(;v) spectrum shows a break in correspondence of 
the rise of this pulse. If we assume for this GRB 
what we found for other GRBs with known afterglow 
emission, this second pulse could be due to an exter- 
nal shock. However the time behavior of Ep is not 
consistent with that expected during the decay phase 
of the fast cooling of a spherical shock. A possible 
interpretation, in the context of an external shock of 
a fireball model, is that during the second pulse we 
are observing the phase in which the Lorentz factor 
of the shocked material is still constant or is decreas- 
ing very slowly. The value of p found for this GRB is 
about 4, that justifies its non detection 17 hrs after 
the primary event. 

We are very grateful to Hara Papathanassiou 
for very useful comments and suggestions on the 
manuscript. This research is supported by the Italian 
Space Agency ASI 
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Table 1 
grbs included in our sample 



GRB 


Position 
error radius 


Offset wrt 
inst. axis 


X-ray ('^^ 
peak flux 


(s) 


7-ray 
peak flux 


T, 
(s) 


TOO 
delay 


Counterparts 


GRB960720 


3' 


11.1° 


0.25 


17 


10 


8 


45*^ 




GRB970111 


3' 


14.6° 


1.4 


60 


56 


43 


16^* 


X ? 


GRB970228 


3' 


13.1° 


1.4 


55 


37 


80 


gh 


X, opt. 


GRB970402 


3' 


8.5° 


0.16 


150 


3.2 


150 




X 


GRB970508 


1.9' 


10.3° 


0.35 


29 


5.6 


15 


5.7^ 


X, opt., radio 


GRB971214 


3.3' 


16.5° 


0.2 


35 


6.8 


35 


6.7^* 


X, opt. 


GRB980329 


3' 


19.2° 


1.3 


68 


51 


58 




X, opt., radio 


GRB980425 


8' 


15.1° 


0.61 


40 


2.4 


31 


gh 


X, opt, radio ? 



Peak fluxes in units of 10 ^ erg cm ^ s 
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Table 2 

Spectra of GRBs and associated X ray afterglow 





c;rl! 


Sc-cl. 


-I'x 


-i\ 


E, 


N// 




() 


1 


GRB960720 


A 


< -0.1 


< -0.1 


> 700 


0.255 










B 


0.39 ± 0.18 


0.39 it 0.18 


> 700 


0.255 










C 


0.49 ± 0.17 


1.76 ± 0.32 


178 it 40 


0.255 










D 


0.76 ± 0.21 


2.21 it 0.40 


28 it 17 


0.255 










E 


1.18 ± 0.31 


> 1.64 


< 16.7 


0.255 










F 


1.9 ± 0.6 


> 2.44 


< 3. 


0.255 






2 


GRB970111 


A 


< -0.37 


0.63 it 0.33 


> 500 


0.46 










B 


< -0.4 


1.18 it 0.17 


> 194 


0.46 










^ 


O fiQ -1- O Qft 

-U.Do it U.oD 


1 OQ -1- n oo 

L.zo ± u.uy 


\ 1 70 


U.4D 










D 


-0.92 ± 0.31 


1.72 ± 0.06 


> 115 


0.46 










E 


-0.76 ± 0.39 


1.73 ± 0.06 


> 118 


0.46 










F 


-0.82 ± 0.27 


1.89 ± 0.05 


> 89 


0.46 










G 


-0.96 ± 0.29 


2.38 ± 0.11 


37 it 12 


0.46 










TT 

xl 


-U.o ( =t U.z / 


o TV 1 nil 


00 1 Q 


U.4D 










I 


-0.19 d= 0.17 


2.9 ± 0.2 


36 ± 9 


0.46 










J 


0.33 ± 0.12 


3.05 ± 0.27 


35 ± 7 


0.46 




> 1.5 


Q 
O 


vjixriy / Uzzo 


A 

A 


o oo _L n OQ 
U.yz it U.Uo 


1 -L n 1 Q 
1.04 it U.lo 


\ 700 


1 










B 


1.4 it 0.1 


2.5 it 0.1 


35 it 18 


1.6 










C 


1.8 ± 0.1 


1.8 ± 0.1 


< 2 


1.6 










D 


1.84 it 0.09 


1.84 it 0.09 


< 2 


1.6 










E 


1.92 it 0.15 


1.92 it 0.15 


< 2 


1.6 










F 


1.5 it 0.4 


> 0.6 


< 2 


1.6 










G 


1.6 it 0.1 


1.4 it 0.3 


< 2 


1.6 


2.06 it 0.24 3.5 it 2.3 


1.33 it 0.12 


4 


GRB970402 


A 


1.38 it 0.03 


1.38 it 0.03 


>700 


< 20 










TD 
ID 


l.oo it U.U4 


1 Qfi J- o rtA 
l.oo it U.U4 


^700 
^ / uu 


1 


i. / it U.O Z.l 


1 J- o OQ 
i.OO it U.Uo 


5 


GRB970508 


A 


0.83 it 0.11 


0.83 it 0.11 


> 700 


0.51 










B 


1.54 it 0.10 


1.54 it 0.10 


>700 


0.51 










C 


1.74 ± 0.14 


1.74 ± 0.14 


>700 


0.51 










D 


1.8 it 0.4 


> 1.4 


< 24 


0.51 


2.1 ± 0.6 10. it 5. 


1.1 it 0.1 


6 


GRB971214 


A 


0.37 ± 0.23 


1.4 it 0.03 


> 700 


0.16 










B 


0.33 ± 0.27 


1.0 it 0.07 


> 224 


0.16 










C 


0.96 ± 0.09 


2.6 ± 0.5 


56 ± 10 


0.16 


1.7 it 0. 2 0.9 it 0.5 


1.2 ± 0.13 


7 


GRB980329 


A 


< 0.16 


1.32 it 0.03 


> 229 


120 it 80 










B 


0. ± 0.33 


1.30 it 0.06 


> 168 


120 it 80 










C 


0.36 ± 0.19 


1.25 ± 0.08 


> 210 


120 it 80 










D 


0.67 it 0.11 


1.29 it 0.02 


> 320 


120 it 80 










E 


0.73 it 0.18 


1.61 it 0.04 


> 213 


250 it 100 










F 


0.78 it 0.32 


1.70 it 0.09 


> 180 


85 it 40 










G 


1.0 it 0.4 


1.7 it 0.2 


> 175 


10. 










H 


1.5 it 0.5 


2.3 di 0.3 


105 it 80. 


10. 


2.4 it 0.4 10 it 0.4 


1.35 it 0.03 


8 


GRB980425 


A 


< 0.3 


1.75 it 0.15 


> 68 


0.396 










B 


0.78 it 0.27 


2.3 it 0.1 


68 it 40 


0.396 










C 


1.15 ± 0.75 


3..S ± 0.7 


.■^1 ± 25 


0.396 










D 


1.0 ± 0.5 


.'is 


-jT 









Rc-r, 



ID) 



^Feroci et al. 1998 , ^Frontera et aL 1998 , ^Nicastro et al. 
Fiume et al. 1999, ^in 't Zand et al. 1998 , ^Pian et al. 1999a 



1998 , "'Piro et al. 1998 , ^Amati et al. 1998 , ''Heise et al. 1998 , ^Dal 



^The photon indices reported in columns 2 and 3 refer to fit with the smoothed broken power-law proposed by Band et al. (1993) 
^Ep values are in keV; Nh values are in units of lO^'^ cm~^ 

°Data referring to the first impulse of the event (see text and paper from Frontera et al. 1998a) 
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Table 3 

grb and x-ray afterglow energetics 



ORB 


Sa- 


S-, 




Sa'/S-, 


S„ /S-, 


S„/Sa 




z 


Rcf.(<^^ 




(2 10 krV) 


( 10 700 krV) 


(2 10 kvV) 


{%) 












GRB960720 


0.08 ± 0.02 


2.6 ± 0.3 




3.0 ± 0.8 












GRB970111 


1.6 ± 0.1 


43 ± 3. 


< 1.1 


3.7 ± 0.2 


< 3 


< 0.69 








GRB970228 


2.2 ± 0.5 


11 ± 1. 


2.16 ± 0.81 


20 ± 0.81 


19.6 ± 0.2 


0.98 ± 0.19 


0.25 


0.695 


(1) 


GRB970402 


0.4 ± 0.04 


8.2 ± 0.9 


0.47 ± 0.15 


4.8 ± 0.1 


5.7 ± 0.7 


1.17 ± 0.13 






GRB970508 


0.33 ± 0.05 


1.8 ± 0.3 


0.47 ± 0.31 


38.9 ± 1.8 


26 ± 12 


0.67 ± 0.30 


~0.078 


0.835 


(2) 


GRB971214 


0.19 ± 0.03 


8.8 ± 0.9 


0.35 ± 0.14 


1.2 ± 0.1 


3.9 ± 0.7 


3.20 ± 0.75 


~0.936 


3.418 


(3) 


GRB980329 


0.7 ± 0.12 


55. ± 5. 


0.72 ± 0.27 


1.0 ± 0.04 


0.40 ± 0.07 


0.39 ± 0.8 


~3.86 (?) 


~5(?) 


(4) 


GRB980425 


0.78 ± 0.02 


2.8. ± 0.5 


<0.46 


28 ± 5 


<16 


<0.59 


< 7.45x10-8 


0.008 


(5) 



^AU fluences are given in units of 10"^ crg/cm^ 
''E62 = afterglow energy in units of 10^^ erg 

"^References for redshift: ^Djorgovski et al., 1999 , "^Mctzgcr ct al., 1997 , ^Kulkarni et al. 1998a , ''Pruchter, 1999 , ^Galama et aL, 
1998b 



TABLE 4 

Derivation of external shock physical parameters 



GRB 


to(s) 


70 


es (early) 




3 arly) 


F,y,rna2; (early) 


eB(late) 


F.,™ar.(late)(^) 


p(b) 


p(b) 




ad rad 




ad rad 


ad 


rad 


(mJy) 




(mJy) 


(t„<t<tc) 


(t>tc) 


GRB970228 


42 42 


110 


0.005 0.0052 


0.61 


0.63 


> 10 


0.001 


7 


2.77 


2.44 


GRB970508 


1000 550 


181±120 


0.046 


0.57 




> 1 


0.037 


~ 2 


2.13 


2.46 


GRB971214 


^IxlC* ~3xl0'* 


180±78 


0.15 




0.65 


~ 0.3 




0.02-0.09 


2.27 


2.60 


GRB980329 


~6.9xl0S 1.5x10^ 


215±80('=) 


0.053 


~1 




-0.35 


0.071 


-0.3 


2.47 


2.80 


GRB980425 




<50 



















^References for F^^jnax(ia-te): Masetti et al. 1999 (GRB970228), Galama et al. 1998a (GRB970508), Dal Fiume et al. 1999 and Ramaprakash 
et al. 1998 (GRB971214), Palazzi et al. 1998 (GRB980329) 
^The typical error on the reported p values is 0.15 
"^Assuming the redshift proposed by Fruchter 1999 
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10 
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40 60 
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Fig. 1. — Light curves of GRBs in our sample, after background subtraction, in the 2-26 keV (upper panels) and 40-700 keV 
(lower panels) energy bands. Horizontal dotted line gives the background reference of each GRB. Vertical dashed lines indicate 
the temporal sections in which we performed the spectral analysis. 
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Fig. 1-continued. — 
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Fig. 1-continued. — 
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Fig. 2. — Spectral energy distributions [vFiy)) of GRBs in the 2-700 keV. For each GRB, in the different panels are rep orted 
the sp ectra in the temporal sections in which we divided the GRB light curve. Solid lines give the best fit of the Band law ( Band 
et al. .993) to the data. Dashed lines show the limit slope in the case of an optically thin synchrotron model (see text). 
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Fig. 2-continued. — 
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Fig. 2-continued. — 
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Fig. 2-continued. — 



F. Prontera et al. 



25 



GRB980329 



1000 



o 

9) 
m 
\ 
N 

B 
o 

> 




100 10 
Energy (keV) 



Fig. 2-continued. — 
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Fig. 3. — Scatter plot of 7-ray photon index (F-y) vs. X-ray photon index (Fx). The dashed hne shows the locus of F-^ = F 
(see text). 
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Fig. 4. Low-ciicrgy power law index Fx vs. peak energy Ep of the v¥(i>) spectrum, for each of the temporal sections of GRBs 
in our sample. The dashed line corresponds to the maximum power law photon index below the energy break, that is consistent 
with an optically thin synchrotron shock model, while the dotted line corresponds to the power law index when cooling of the 
high energy electrons is taken into account (see text). The dashed-dotted line correponds to the limit photon index of the Inverse 
Compton (see text). 
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Fig. 5. — Ratio between the 2-10 keV fluence derived from the GRB Mght curve (GRB fluence) and that derived from the 

late afterglow observation (afterglow fluence) vs. the temporal index 5 of the late afterglow fading law. The afterglow fluence is 

integrated from the end (T) of the GRB light curve, as estimated by our data (see Table 1) to l(f s (see text). 

Top: The GRB fluence is integrated over the time interval starting from 1% of the GRB time duration to its end. 

Bottom: The GRB fluence is integrated over the time interval starting from 63% of the GRB time duration and the end of the 

event. 

Superposed to the data are the expected ratios assuming the GRB fluence has the same origin as the late afterglow emission. 
Continuous line: best fit curve; dashed line: expected curve if there is no spectral softening from the prompt to the late afterglow 
emission. The best fit curve in the top case would require a hardening of the spectra that is in contrast with the observations. 
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Adiabatic Shock 
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Fig. 6. — Theoretical afterglow spectra at various times from the afterglow onset, expected on the basis of the synchrotron shock 
model by Sari, Piran and Narayan (1998). Top: a fully adiabatic shock; bottom: a fully radiative shock. Triangles correspond to 
the break energy Em, while the squares correspond to the cooling break Ec (see text for definitions). The dotted vertical lines 
limit our energy passband. The parameters assumed are: energy of the shock E = 10^^ erg, index of the power-law distribution of 
the electrons p = 2.3, particle density of the interstellar medium n = 1 cm~^, fireball distance D = 10^* cm, initial Lorentz factor 
of the ejecta 70 =100. The other parameters (see text) ee and es depend on the type of cooling of the shock. We have assumed 
£b = 0.01 and ee = 0.1 for an adiabatic shock, es = 0.1 and ee = 0.8 for a radiative shock. 
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Fig. 7. — Top panel: High-energy power-law photon index during the GRB tail vs. index of the afterglow fading power-law. 
Bottom panel: photon index Fa of the power law spectrum of the X-ray afterglow vs. index of the afterglow fading power-law. 
Dotted line: expected relationship for a fully adiabatic shock if the afterglow observation is performed at time t>tc; DoUed-dashed 
line: expected relationship for a fully adiabatic shock if the afterglow observation is performed at time t, with tm <t<tc; Dashed 
line: expected relationship for a fully radiative shock. See also text. 
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Fig. 8. — Time behaviour of the measured peak energy Ep of the uF{i/) spectrum for GRBs in our sample. Vertical dotted line- 
time at which early afterglow is expected to have already started (see Sect. 5.2). Ep values that are expected to be coincident with 
Fim are shown as full circles, while those coincident with Ec are shown with open triangles. The expected time behaviour of Ep in 
the case of the synchrotron shock model by Sari, Piran and Narayan (1998) is also shown for two solutions of the shock: adiabatic 
(continuum lines) and radiative (dashed lines) cooling. Line with higher slope corresponds to fast cooling and that with lower 
slope corresponds to slow cooling. See text for details. 
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Fig. 8-continued. — 



